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Abstract 

We investigate the violation of time reversal invariance in the decay of the free neutron in 
the framework of the Minimal Supersymmetric Standard Model (MSSM). The coefficient of 
the triple product of the neutron spin and the momenta of electron and neutrino, the so called 
D parameter, is computed at one loop order including all diagrams. We find that D is mainly 
sensitive to the trilinear A coupling in the squark sector and to the phase of the coefficient \i 
which mixes the two Higgs superfields. The maximal MSSM contribution using parameters 
still allowed by experiment is however at D « 10~ 7 , while QED final state interactions give 
a value of _Df s ; = —2.3 • 10~ 5 . Explicit expressions for all relevant diagrams are given in an 
Appendix. 



1 Introduction 



The existence of complex parameters in the Lagrangian can lead to a violation of time reversal 
symmetry. In the standard model the phase of the CKM matrix is the only phase which cannot 
be eliminated by a field redefinition. In the MSSM there are additional parameters which in 
general cannot be made real: 

• the coefficient of the term bilinear in the Higgs superfields: \i 

• two of the three gaugino masses rhi i — 1 , 2 

• terms mixing left- and right-handed sfermions Af. 

The violation of (naive) T symmetry can be tested using an observable which is odd under 
applying the T symmetry operator. For example an odd combination of spins and momenta 
fulfills this condition. In this paper we investigate the triple product 

— • (Pe x Vv) (1) 

of the spin of the neutron and the momenta of electron and electron antineutrino in the decay 
of free neutrons. The coefficient of this expression in the decay distribution 



d-Egdcos^ep t a n E e E p 

where G e is the tree level expression, is called the D parameter. It offers a distinct possibility 
to search for T symmetry violation in neutron decay, besides the electric dipole moment of the 
neutron. An experimental effort is currently underway at the ILL Grenoble to improve the 
measurement of, or bound on, D. A complete calculation of D in the MSSM therefore seems 
timely. To the best of our knowledge, only the gluino loop contribution to D has previously 
been calculated PQ. 

In this paper we follow the conventions of Rosiek [2] ; that reference also contains expressions 
for all relevant Feynman rules. A brief summary of the (somewhat unusual) notation is given 
in Appendix A. 



2 D— Parameter 
2.1 Standard Model 

As remarked earlier, the only T— violating parameter in the Standard Model (SM) is the 
phase in the Kobayashi-Maskawa matrix. It can lead to violation of T (or CP) symmetry 
only in processes involving all three generations of quarks. Therefore it can contribute to the 
D— parameter only starting at the two-loop level. As a result, the truly T symmetry violating 
contribution to D is very small in the SM 

D SM < 10- 12 . (3) 

Experimentally a complete time reversal, which would consist of motion reversal and ex- 
change of the initial and final states, is unfortunately not possible in neutron decay. Instead, 
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D is odd under so-called naive time reversal, where only the directions of all spin and momen- 
tum three-vectors are reversed, without exchanging initial and final state j3]. While genuine 
T invariance can only be violated if some parameters in the fundamental Lagrangian contain 
nontrivial complex phases, naive time reversal invariance can be violated whenever the rele- 
vant matrix element has a nonvanishing imaginary part. This difference is significant, since 
an imaginary part in the matrix element, a so-called absorptive phase, can also originate from 
loop corrections which respect genuine T invariance. In the present case these are due to 
QED final state interactions between the proton and the electron. Note that a loop diagram 
gives an absorptive phase only if the particles in the loop can be on-shell; this leads to an 
additional phase space suppression factor of order E e /m n , where m n is the neutron mass and 
E e the energy of the electron in the neutron rest frame. The total contribution from final state 
interactions is therefore quite small Jj]. 



the bound is saturated at the kinematic maximum of E e . "New physics" contributions to 
D that are much smaller than this value will be very difficult to extract even for arbitrarily 
small experimental error, since the prediction (£Q) has some theoretical uncertainties, e.g. due 
to higher order corrections and proton form factor effects. Finally, the current experimental 
sensitivity jB] is still well below the prediction (J3J, 



However, efforts are underway to improve the sensitivity by nearly an order of magnitude [4. . 



Let us now turn to the calculation of the D parameter in the MSSM. We have extended the 
analysis of ref. 1J by including all possible diagrams at one loop order. Four different types of 
diagrams can contribute to neutron decay: 

• vertex correction at the W— quark vertex; 

• vertex correction at the W— lepton vertex; 

• vertex corrections where the exchanged W boson is replaced by a charged Higgs boson; 

• box diagrams. 

The corrections to the W— lepton vertex give a contribution to D that is suppressed by a 
factor — ~ 5 ■ 10~ 4 , so these diagrams can safely be neglected compared to the corrections 
to the W— quark vertex. The diagrams with Higgs boson exchange do not contribute at all to 
the D— parameter, since they do not contain sufficiently many 7 matrices to give rise to a spin 
correlation. We therefore only need to consider corrections to the W— quark vertex as well as 
box diagrams. 

Since we are computing a contribution to the neutron decay distribution which has non- 
trivial dependence on final state momenta, see eq.(J2J), we cannot completely neglect external 
momenta when evaluating the loop integrals, even though these momenta are much smaller 
than the masses of the superparticles in the loop. However, after introducing Feynman pa- 
rameters and shifting the loop integration variable k in such a way that the terms linear in 



D M \ < 2.3 ■ 10 



;-5. 



(4) 




(5) 



2.2 MSSM 
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k are eliminated from the denominator, all terms of order m n or m e can be neglected in the 
denominator; in other words, external momenta can be ignored in the loop integrals after the 
shift of the loop momentum. The three- and four-point functions that appear in our calcu- 
lation can therefore easily be reduced to combinations of two-point functions, as described in 
Appendix C. 

The coefficients in front of the loop integrals contain three different kinds of suppression 
factors. The Dirac algebra can introduce factors of the nucleon mass, rather than the mass of a 
fermionic superparticle. Moreover, certain chargino and neutralino couplings contain Yukawa 
couplings to first generation fermions. Finally, a term may require mixing between SU(2) 
doublet and singlet first generation sfermions, which is again proportional to a first generation 
Yukawa coupling. Numerically these three suppression factors are of comparable size, so that 
a simple counting of these factors is sufficient to isolate the leading terms. 

Let us illustrate these remarks by analyzing the k° — u — d loop correction to the Wud 
vertex, see Fig. 10 in Appendix El The corresponding contribution to the D parameter is given 
in eq. (|B.5|) . We first note that the three-point function C^ u is 0(1), whereas the functions 
Co, Cd2 and Cd3, which are defined in Appendix C, are C(l/mg USY ), where msusY ~ 0.1 to 
1 TeV is a typical superparticle mass scale. Let r = m„/m SU sY, and Y u and Y d be the u and d 
quark Yukawa couplings, respectively. From eqs. (|B.6|) for the relevant couplings one can then 
derive the following behavior for the various terms listed in eq. (|B.5|) . which we label here by 
the relevant product of couplings: 

A 1 D 1 E 1 ^0{Y u Y d ) ] A l C 1 E l ^O(rY u ); B 1 D 1 E 1 ~ 0(rY d ); B 1 C 1 E 1 ~ 0(r 2 ). (6) 

We note that each term has two suppression factors. This is true for all other loop corrections 
as well, as can be seen from the results given in Appendix B. Terms with an even larger 
number of suppression factors have been omitted. 

It should be noted that the terms containing powers of the suppression factor r have been 
obtained by replacing the kinematic u and d quark mass by the mass of the proton and neutron, 
respectively; this simple approximation goes under the name of "naive dimensional analysis" 
[7]. It seems reasonable to take some sort of long-distance quark mass here, although the 
use of a constituent quark mass ~ m n /3 could also be defended. Since the first generation 
Yukawa couplings are* (9(10~ 4 ) one might think that Y ud give a much more severe suppression 
than r. However, this need not be the case. Only the imaginary parts of the products of 
couplings are relevant. Terms with explicit factors of Yukawa couplings can acquire nontrivial 
phases either from gaugino-higgsino mixing in the chargino-neutralino sector, or from mixing 
between SU(2) doublet and singlet sfermions, whereas (in the absence of sflavor mixing) terms 
without Yukawa couplings are generally not sensitive to phases in the sfermion sector. Note 
also that Y d oc tan/3 for tan/3 ^> 1. These considerations imply that contributions that are 
suppressed by Yukawa couplings are in general not smaller than those that are suppressed 
by powers of r. Finally, the box diagrams receive an additional (modest) suppression factor 
of order (mw / '"^susy) 2 • The same factor should also be multiplied to all terms that do not 
contain factors of r, while terms with only one power of r receive additional suppression of 
order m w / m SU sY • The reason is that heavy superparticles must decouple quadratically. Such 
factors result from mixing in the k sector, and/or from mixing between SU(2) doublet and 
singlet squarks. [The latter is also suppressed by the relevant Yukawa coupling, but this is 
already included in eq.©-] However, for reasonable superparticle masses these additional 
suppression factors are much less important than the factors listed in eq.©. 

*These are short-distance couplings, to be taken at a momentum scale of order TOsusy- 
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2.3 Restrictions on the parameter space 

In order to make a useful analysis of the D parameter it is important to know which parts 
of the parameter space of the MSSM are still experimentally allowed. Especially relevant are 
experiments which test T symmetry violation in other observables, in particular the electric 
dipole moments (EDMs) of electron and neutron, d e and d n . Current experimental bounds on 
these quantities [6] impose strong constraints on parameter space. For superparticle masses 
of order (a few) hundred GeV, many combinations of phases are excluded, although in some 
cases there is a possibility that one can have small EDMs while retaining large phases jSJ. 
For our analysis we have used the formulae for the EDMs given in [Hj, and checked for each 
parameter point that it does not violate the experimental limits 

|4| < 0.63 ■ l(T 25 e cm (CL=90%) 
-0.005 ■ 10~ 26 e cm < 4 < 0.143 ■ 10~ 26 e cm (CL=68%) (7) 

Of course, limits on superparticle masses from null results of experimental searches for these 
particles at high energy colliders also have to be respected jH]. 

3 Numerical Analysis 
3.1 Choice of Parameters 

We wish to find the maximal supersymmetric contribution to D in the framework of the 
R— parity conserving MSSM. As well known, the parameter space of the general MSSM is 
vast, so some simplifying assumptions are necessary. In our analysis we have assumed that no 
flavor mixing exists in the sfermion sector. This means that the ffH and ffg vertices are 
diagonal in flavor space. One can then easily see from the diagrams given in Appendix B that 
CKM mixing between quarks is not relevant, i.e. all sfermions in the loop must be of the first 
generation. On the other hand, our discussion of eq. © showed that it is of crucial importance 
to include L — R mixing even for first generation squarks. 

Ignoring all flavor mixing between sfermions may sound like a rather severe restriction on 
the parameter space. However, mixing between first generation sfermions and those of the 
second or third generation is strongly constrained by experimental limits on various FCNC 
processes [5]. We have checked that including flavor off-diagonal LR and RL entries of the 
experimentally allowed magnitude in the squark mass matrices does not increase the maximal 
contribution to D once the constraints (J7J) on the EDMs have been imposed. Generally 
speaking, flavor mixing should not be very important here, since all external fermions belong 
to the same (first) generation. 

In the absence of flavor mixing between sfermions our results are independent of the soft 
breaking parameters describing the second and third generation sfermion mass matrices. The 
same is true for the soft breaking parameters of the tree-level Higgs potential. We have 
initially chosen a rather small value for the ratio of vacuum expectation values (vevs) tan/3, 

tan/5 = 3. (8) 

The reason is that supersymmetric contributions to the EDMs increase with tan f3 jH] • Choos- 
ing a small value for this parameter therefore minimizes the impact of the experimental con- 
straints (J7j) on the allowed values of the remaining parameters. 
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Supersymmetry is a decoupling theory, which means that in the limit where the masses of 
the super symmetric particles go to infinity, the predictions for all observables approach their 
SM values. We therefore consider a relatively light spectrum of superparticles, described by 



the following values of the relevant soft breaking parameters: 

m 2 L = 35 ■ 10 3 GeV 2 m\ = 50 • 10 3 GeV 2 (9a) 

m\ = 150 • 10 3 GeV 2 = m 2 D = 200 • 10 3 GeV 2 (9b) 

\n\ = 450 GeV (9c) 

|m x | = 200 GeV |m 2 | = 400 GeV \m 3 \ = 800 GeV (9d) 



Here and are the soft breaking masses for SU(2) doublet and singlet sleptons, vhq 
is the soft breaking mass for SU (2) doublet squarks, and mjj and mjj are the soft breaking 
masses for SU(2) singlet squarks with charge 2/3 and —1/3, respectively; recall that we only 
need to specify these masses for the first generation. As mentioned earlier, /i is the coefficient 
of the term coupling the two Higgs superfields in the superpotential. Finally, rrii are the soft 
breaking gaugino masses; the ratios of these masses are similar to that expected in Grand 
Unified models with universal gaugino mass at the unification scale. The choices © lead to a 
superparticle spectrum that respects all experimental limits from searches for superparticles, 
and allows large CP-violating phases through the cancellation mechanism jH]. We will later 
comment on the effect of lowering the overall SUSY mass scale from the choice of eqs.fjHJ). 

The remaining parameters are sampled randomly, taking a flat distribution within specified 
limits. The trilinear couplings^ 

|4J = |4|, |4| e [0,0.1] GeV (10) 

were chosen so that internal cancellations in the EDMs are possible but the sfermions do not 
acquire vevs, which would be the case for too large values. For simplicity we took the same 
value for and A u ; A\ is allowed to differ, in order to facilitate independent cancellations 
in the supersymmetric contributions to d n and d e . Finally the phases of the gaugino masses, 
[A and of the trilinear couplings were varied independently (except A d = A u ) over the entire 
possible range, 

0m 2 , <pA u = <frA d , <\>A V M G [0,27r[. (11) 

Note that by an appropriate field redefinition the phase of the gluino mass can always be set 
to zero without loss of generality. 

3.2 Numerical Results 

In this section we display the results of our numerical analysis. Note that only about 20,000 
out of 10 10 tested sets of parameters satisfied the constraints (JJJ). This illustrates that the 
EDMs do indeed severely constrain the allowed combinations of nontrivial complex phases in 
the MSSM Lagrangian. 

In Fig. [T] we plot the dependence of the supersymmetric contribution -Dsusy to the D 
parameter on the phase </> M . It is easy to see that -Dsusy depends strongly on this phase. 
There is however also a large variability at fixed M , which shows that -Dsusy also depends 
significantly on the other parameters. 

^Recall that we are using the convention of ref.|2], where the ordinary Yukawa couplings are explicitly 
included in the A— parameters. In case of first generation quarks these couplings are roughly of order 10~ 4 ; 
this explains the small values of the A— parameters in ea. (|10|l . 
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Figure 1: Supersymmetric contribution to the D parameter depending on the phase 

Before investigating this closer it is useful to identify the diagram which gives the leading 
supersymmetric contribution to D. Inspection of the contributing diagrams in Appendix B 
shows that only one of them, the quark vertex correction shown in Fig. 11, involves the 
strong interactions. One might have thought that the presence of the gluino in this diagram, 
which is significantly heavier than the electroweak gauginos, would partially compensate this 
enhancement. However, this would only be true if the gluino was significantly heavier than 
the squarks in the loop; such an ordering of masses is not allowed in the MSSM, since it would 
lead to tachyonic squark masses at energy scales just above the gluino mass [TU] . Indeed 
we find numerically that the gluino vertex diagram, which is the only diagram considered in 
PQ, gives the leading contribution to Dsusy- All other diagrams are suppressed by at least 
one order of magnitude. This is in spite of the fact that the only phases contributing to the 
gluino loop diagram come from squark mixing, whereas the electroweak loop corrections are 
also sensitive to phases from electroweak gaugino-higgsino mixing."'" On the other hand, the 
relative importance of the other diagrams is increased by the fact that they all add destructively 
to the gluino loop diagram, i.e. tend to reduce | -Dsusy I- As a result, the pure gluino loop 
contribution to -Dsusy > shown in Fig. is somewhat larger in absolute size than the total 
one-loop contribution shown in Fig. ^ 

The contributions of the other five diagrams are shown in Fig. El For better readability 
the y— axis was scaled down by a factor of 10. The numbers in brackets in the label of the y 
axis denote the diagram whose contribution is shown in the corresponding plot. 

We see that the first of the three electroweak vertex corrections, with a neutralino and two 
squarks in the loop, gives significantly larger contributions than diagrams 3 and 4, which have 
one squark, one chargino and one neutralino in the loop. This can be understood from the 

^Partly for this reason, the chargino and gluino loop contributions to d n can be of similar size. Note, 
however, that in case of the D parameter the chargino and gluino loop contributions have quite different 
structure, i.e. there is no chargino diagram with two squark propagators. 
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Figure 2: The contribution from the leading gluino diagram to -Dsusy- 

observation that the Wkfk® vertex only couples two SU(2) gauginos (winos) or two higgsinos 
together; this suppresses possible contributions involving the phase fa associated with the 
U(l)y gaugino (bino). Moreover, the coupling structure in diagrams 3 and 4 is such that some 
gaugino-higgsino mixing is required, whereas diagram 1 gets finite contributions even without 
this mixing. This suppresses the contributions of diagrams 3 and 4 by a factor 0(mw/msnsY) 
relative to that of diagram 1. Finally, for the given choice of sparticle masses, the box diagrams 
5 and 6 give contributions which are about two orders of magnitude smaller than that of 
diagrams 3 and 4. We remarked earlier that the contribution from box diagrams are suppressed 
by (^wV^susy) 2 — 1/20 for our set of parameters. Moreover, the additional integration 
over Feynman parameters required in the D— functions appearing in the box contributions 
gives another suppression factor ~ 1/5 compared to the C— functions appearing in the vertex 
corrections. 

Let us now return to the question which other parameters influence the size of .Dsusy- As 
the gluino diagram is independent of the phases <f> mi and (p m2 these cannot play significant 
roles. If we now restrict M to a small interval around it (0 m = 7r ± 0.0l7r), in accordance with 
the limits on the EDMs, it becomes clear that A u is the second parameter which determines 
-Dsusy for given sparticle masses. \A U \ largely determines the amount of L — R mixing between 
up-type squarks, since in this case the contribution oc is suppressed by a factor cot j3\ many 
of the terms oc Y u originate from this mixing. Since \x is almost real, | -Dsusy I becomes very 
small as \A^\ —>■ 0, as shown in Fig. 4. Moreover, for (almost) real fj,, the phase of A u (= A& for 
our choice of parameters) determines the CP-violating phases in the squark mixing matrices. 
This leads us to expect that | -Dsusy I will be maximal if 4>a u is ~ | or ~ This is confirmed 
by Fig. El (| -Dsusy I can be small even for these choices of (f>A u since \A U \ is still varied in Fig. 03 
of course, the value of 4>a u becomes irrelevant as \A U \ — > 0.) 

We have also investigated the influence of tan/3 on .Dsusy- To that end the parameter 
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Figure 3: Electroweak contributions to -Dsusy, The numbers in parentheses correspond to 
the labeling of diagrams and their contributions in Appendix B; e.g. -Dsusy(I) refers to the 
diagram shown in Fig. 10, whose contribution is given by eq. ()B.5j) . 
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was chosen arbitrarily from the set of points that are allowed for tan/3 = 3. The tan/3 
dependence of -Dsusy is shown in Fig. |H1 The almost linear increase results from the increase 
of the d— quark Yukawa coupling, and hence of di — d^ mixing, which is proportional to 
1/ cos/3 ~ tan/3 for tan 2 /3 3> 1. Here we have neglected the restrictions from the EDMs; for 
the set of parameters described by eqs.Q, (fTUj) and (fT2J). these impose the bound tan/3 < 3.5. 
On the other hand, since several contributions to d n and d e grow oc tan/3, cancellations 
can also work at large tan /3, if some of the other parameters, e.g. the phases, are changed 
slightly, without significantly modifying the prediction for .Dsusy- However, since the separate 
contributions to d e and d n grow with increasing tan /3, increasingly precise cancellations become 
necessary to satisfy the experimental constraints (J7J). 

Both of these observations are confirmed by Fig. 7, which shows Dsusy for the same 
parameters as in Fig. 1, except that now tan/? = 30. We see that the overall scale of Dsusy 
is one order of magnitude larger than in Fig. 1, as expected from the linear growth shown 
in Fig. 6. The small number of surviving points illustrates the difficulty of getting both \d e \ 
and \d n \ sufficiently small through delicate cancellations. In particular, the d e constraint now 
excludes some region of M altogether. § Recall, however, that the (s)leptonic corrections to 
D^susy are suppressed by a factor m e /m p , and are thus negligible. We can therefore vary 
the sleptonic soft breaking masses mi, ttir without significantly changing the prediction for 
D^susy- in this case the entire range of M becomes allowed again. We therefore conclude that 

§ Recall that </y measures the relative phase between /j, and the gluino mass m^. For the choice © of 
dimensional parameters the d e constraint limits the relative phase between /i and the SU (2) gaugino mass mi 
to narrow bands around and 7r even for small tan/3, but this is not visible after scanning over </> TO2 . 
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Figure 7: -Dsusy vs. for tan/3 = 30. 



the maximal allowed | -Dsusy | increases essentially linearly with tan/3. However, requiring the 
bottom Yukawa coupling to be less than that of the top quark, or at least to be sufficiently 
small to not have a Landau pole below the scale of Grand Unification, leads to the upper 
bound tan/3 < 60. If the dimensionful parameters in the squark, gaugino and higgsino sectors 
are as in eqs.(9), the maximal value of |-Dsusy| is therefore around 5 • 10~ 8 . 

Our choice J0J) for the relevant soft breaking masses means that first generation squarks as 
well as most charginos and neutralinos have masses around 400 to 500 GeV. Direct searches 
for sparticles allow us to lower these masses by a factor of two to three. We expect that for 
fixed phases D scales quadratically with the overall superparticle mass scale, D oc l/mg USY - 

This is borne out by Fig. 8, which shows -Dsusy for tan/3 = 3; the other parameters are 
as in eqs.(jnj and ()12)) . except that all quantities with mass dimension are multiplied with the 
dimensionless factor c, which is varied between 0.4 and 2. Similar to the case of Fig. 6, values 
of c significantly different from unity are disallowed by the constraints on d e and/or d n , but 
again this can be fixed by small variations of the phases. Since -Dsusy and the dipole moments 
show the same c~ 2 dependence on c, increasingly delicate cancellations are required to satisfy 
the constraints ((7J) as c is reduced. The lower bound c > 0.44 is in our case set by the lower 
bound on the mass of the lightest selectron, m g > 95 GeV, from LEP searches 

By simultaneously reducing the overall SUSY mass scale and increasing tan f3 one might 
therefore in principle be able to reach values of |-Dsusy| slightly above 10~ 7 . However, the 
finetuning required to satisfy the limits on both d e and d n then becomes very severe indeed; less 
than one parameter set in 10 9 will survive.^ Moreover, one would have to choose soft breaking 
masses for f sleptons and b squarks that are much larger than those for the corresponding 
first generation sfermions. Otherwise L — R mixing, which grows oc |/i|tan/3, would make 
the lighter f and b mass eigenstates much too light, or even tachyonic. Recall that even 
I-DsusyI — 10~ 7 is still two orders of magnitude below the contribution (j^J) from final state 

"Since d e and d n have to be finetuned independently, the overall severity of finetuning scales like c -4 - tan 2 (3. 
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Figure 8: -Dsusy for tan/3 = 3, where the other parameters are as in eqs.(jOJ) and (fT2j). except 
that all mass parameters have been multiplied with the dimensionless scaling factor c. 

interactions. 

Finally, we have checked for a few cases that changing the ratios of soft breaking parameters 
from the choice of eqs.(P|) does not increase the maximal allowed value of .Dsusy significantly. 
This is not surprising, since the overall scale of -Dsusy is set by the heaviest superparticle 
that occurs in a given loop diagram, whereas the lower bound on the overall mass scale is 
essentially set by the lightest (charged) superparticle. |D ) susy| will therefore be maximal if the 
parameters are chosen such that the mass splitting between superparticles is relatively small, 
which is true for the parameters of eqs.©. 

4 Conclusions and Outlook 

In this paper we have analyzed T symmetry violation in the beta decay of free neutrons via 
the D parameter. We have extended the analysis of pQ by computing all diagrams that occur 
at one loop order. We have performed a full scan of the allowed phases and the magnitude 
of the A parameters describing mixing between SU(2) singlet and doublet squarks, subject 
only to the experimental constraints on the electric dipole moments of the electron and, in 
particular, the neutron. 

We find that the gluino loop correction to the Wud vertex indeed gives the leading super- 
symmetric contribution to D. The phase of fi, which has been neglected in ref.pQ, crucially 
influences the size of _D SU sy- Moreover, allowing for cancellations between various contribu- 
tions to the EDMs permits larger values for the relevant phases, which increases | -Dsusy |- We 
nevertheless find that the maximal contribution to \D\ from the R— parity conserving MSSM is 
at least two orders of magnitude smaller than the contribution (j3J) from electromagnetic final 
state interactions. This means that even greatly improved experimental upper bounds on \D\ 
will not lead to new constraints on this model. On the other hand, a measurement of D which 
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differs significantly from the prediction (0J) would rule out the R— parity conserving MSSM 
along with the SM. Larger contributions might be possible in the R— parity violating version 
of the MSSM. Note that R— parity violation only through trilinear terms in the superpotential 
does not contribute to the EDMs at one-loop level [H], whereas e.g. baryon-number violating 
(A") couplings can contribute to D at the one-loop level. 

Analogous triple products can also be defined for decays involving heavier particles than 
the neutron. Decays of heavier baryons, e.g. A or can be treated using the expressions 
given in Appendix B, since here the external momenta are still much smaller than mgusY- 
In these cases the supersymmetric contributions are expected to be larger by several orders 
of magnitude than in case of neutron decay, since Yd would be replaced by Y s or even Yj,; in 
addition, the constraints on CP-violating phases of soft breaking parameters in the second and 
third generation are much weaker than for the first generation. Experimental measurements 
will probably be difficult in these cases, however. 

Even larger supersymmetric contributions can be expected for the analogous asymmetry 
in top quark decay. Since the top Yukawa coupling is 0(1), electroweak corrections might 
well be comparable to SUSY QCD corrections in this case. Moreover, the mass ratio 
m t/ m svsY is also 0(1). This again increases the level of the expected corrections; it also means, 
however, that the expansion for small external momenta used in our calculation is no longer 
adequate. Note also that the spin of the top quark cannot be measured directly; nevertheless 
the D— parameter in top decays does contribute to measurable T— odd asymmetries [12... The 
final state interactions for such decays have already been found to have approximately the 
same size as for neutron decay jTHj. A full calculation of supersymmetric contributions to CP- 
violation in top decay, including electroweak contributions, might therefore prove rewarding. 
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For the benefit of the reader we have summarized in this appendix the conventions used in 
this paper. A complete description including all expressions for the Feynman rules can be 
found in 

There are two charginos K,f,i = 1,2 whose mass matrix is diagonalized by two unitary 
matrices Z + and Z_: 



This equation does not specify the two matrices Z + and Z_ uniquely. This can be used to 
choose both masses positive and sorted in ascending order. 

Similarly, the neutralinos are denoted by i — 1 . . . 4 and the neutralino mass matrix is 



A Notation 
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diagonalized by a unitary matrix such that 



/ mi 





—evi 
2c w 


ev2 
2c w 


\ 





m 2 


evi 
2sw 


— ev2 
2s\v 




— evi 
2c w 


evi 
2sw 





-fl 




\ ev 2 


~ev2 
2s\v 


-n 





/ 



Z N = diag(m K o,...,m K o) 



(A.2) 



As noted in the text, the mass parameter of the 577(3) gauginos can be taken as real. 
Therefore the eight gluinos all have a mass of m.3. 

Finally, the sfermion mass matrix can be written compactly for all four different types of 
sfermions as 



M) 



m 2 j T + m 2 f + 



■Qf*w) 



■Is 



Hcrc. 



f 

rrif (k/j, + Y 
cot j3 for up-type squarks and k 



2 

m ? 



nif \ Kfi + 
+ m 2 f + Q f - 



Y 



(A.3) 



" 1 



"fa ' ■ f 

tan j3 for down-type squarks and charged 



sleptons. m; and denote the mass parameters and Af the coefficient of the trilinear 
terms from the soft SUSY-breaking terms in the Lagrangian, and Y is the respective Yukawa 
coupling. Qf is the electromagnetic charge and Tj? the quantum number of the third component 
of the isospin operator. As there is no SU(2) singlet sneutrino only the upper left element in 
the matrix of eq. IA. 31 occurs for the sneutrinos. 

These mass matrices can be diagonalized with a unitary matrix each, yielding 



z.Miz: 



diag (m~ 
diag (m 



vv ■ 

2 

U\ 1 



Z[MjZ L = diag (m 



m 



11 V, J 



Z ] D M\Z D 



diag 



(A.4) 



B Feynman diagrams 

Altogether six Feynman diagrams were computed. The vertex corrections at the W-lepton 
vertex are, as already mentioned, suppressed by a factor of ^ ~ 5 • 10~ 4 and can therefore be 
neglected. We first give the tree level expression for the differential decay distribution, since 
it enters the definition of the D— parameter: 

D _ dTj ( dr trec f pe x p D \\ 

dE e dcos9 eP 1 Vd£ e dcos# eP ' Sn ' \ Ee E D ) ) 1 ' ' 

All necessary traces have been computed with the help of FORM |14| : in many cases the 
results have been checked using manual calculations. 



Tree level result 



dr ~ - EtmnUV (B . 2) 



&E e & cos 9 eD 87r 3 m^ (m n - E e - E p 
with 



U = C IJ (B.3) 



v^sintf 



V 2 sm Uw 



15 



p L=l e K=l 

Figure 9: Tree level diagram describing neutron decay. 

Here e is the QED coupling constant, 9 W the weak mixing angle, and C the quark flavor 
mixing matrix. In our numerical calculation we have ignored all flavor mixing. 

Diagram 1: Neutralino-w-c/ vertex correction 




Ve e 

Figure 10: Vertex correction with neutralino-w-d loop. 
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with 



23m {AxDxEx) m p )C^(0, 0, m^, m^m n o)8^ 



(B.5) 



+ f43fm(Ai(7i£?i) +43m(B l D l E 1 )jm- K om z n C D2 (0,0,m- K o,rnu. 



+ 43m(B 1 C 1 E 1 )m n C fl 2(0,0,m^o,m ii ,mi) - C D3 (0, 0, m s o, m Si , 



+ . . . 



2^6 



' e -7(7+3)1 ry\l V I ryAl ryli 

— — — — Zj v /j n - Y u ^ N ^u 

6 cos ta 



7 \v 



(B.6a) 



Ci = - 



e 



e 



Z# sin0 w + Z^* cos^J - YlZfZp^ (B.6b) 

(B.6c) 



z5(^sin^-^cos%)+F/4 J+3) ^ 
sin t»V cos "w ^ ^ 

-V2e y(J+3)j rv\U 1 yj 7 Jj y3l 



Ml 

e 



v^sinfl 



D U ^ 



(B.6d) 
(B.6e) 



The dots (...) denote additional terms in the decay distribution that do not contribute to D. 
Eqs. (1B.5J) and ()B.6j) . as well as all subsequent expressions, are given in the notation of Rosiek 



Diagram 2: Gluino-tt-<i vertex correction. 



dr 2 



d-ELd cos d-ELdcos^, 



dr tree _ f p e xp p \ 4 1 <Rt(UV r < 



EpEry J 34mr,U 2 V 2 16tt 2 



-23m (A 2 D 2 E 2 ) m p )C^(0, 0, m Si , m d ~ m^)^ 



+ (4^m(A 2 C 2 E 2 ) +4arm(B 2 D 2 S 2 ))mg ! m^C /J2 (0,0,m £ 5 ! ,m fii ,m d - 



(B.7) 



+ 43m (B 2 C 2 E 2 ) m s n C D2 (0,0 mi)- ^(0,0 



+ . . . 
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with 



A 2 
B 2 



g 3 V2zi; +3)t 

-Q3V2ZH 



C 2 = -g 3 V2Z J D j 

Jj* ryli* (jU 



77, _ 7 JJ* rv 

& 2 — Zj d Zj v 



(B.8a) 
(B.8b) 
(B.8c) 
(B.8d) 
(B.8e) 



Diagram 3: Neutralino Chargino— d vertex correction 



jj=i 




1=1 



Figure 12: Vertex correction with neutralino-chargino— d loop. 
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dr. 



dr 



dE e d cos dev dE e d cos 9 



tree -> I Pe x Pv 

■ S. 



mt(uv 2 ) 



E e E p J 2m n U 2 V 2 16tt 2 



-3m (A3C3F3) m p J (7^(0, 0, m- K o , m~+ , m^ v 

+ 2( 3m m p m-+m«oC o (0, 0, m;,m k o,m k +) 

+ 3m (B 3 D 3 E 3 ) m 3 n C D3 (0, 0, m;,m s o,m s+ ) 



(B.9) 



+ 



(jm (A3D3F3) m-+ + Jm (B 3 D 3 F 3 ) m k ^j m 2 t C D2 (0, 0, mj.,m«o,m-+)j 



+ 



with 



a 3 =y tt J z# i *z?'c 



B, = 



S1116 1 



_r?Mi* glj* _|_ r?(M+3)i* ^2j* \ qIM 



sm.9 w V V2 



£"3 



1^-0 ^21* _j_ _ y 2j r?3l* 

sin9 w \ ~ N ^2 



1 



—7= . 2 2~ Z D Sil1 ° W ~ Z N COS ^ 

V 2 smt^ costly \«J 

7 (J+3)i 7 ll* . V J 7 Ji 7 3l* \ ?JM 
N ' * d D ^ N ]° 



3 COS 9 W 



(B.lOa) 
(B.lOb) 

(B.lOc) 

(B.lOd) 



y/4 J+3)i Z$U JM (B.10e) 



(B.lOf) 



Diagram 4: Chargino-Neutralino-M vertex correction 



dr 4 



dr 



tree 



d_E e d cos 9 eP dE e d cos 9 eP 



Pe X p p 



mt(UV 2 ) 



E e E p J 2m n U 2 V 2 16tt 2 
-3m (A4C4F4) m^C^O, 0, m- K -, m k o, m fii )5 /u/ 

+ 2 (jm m p m k om k -C (0, 0, m fii , m--, m^o) 

+ 3m (B A D A E A ) mlC D3 (0, 0, m Si , m~-, m s o) 



(B.ll) 



+ (3m (A4D4E4) m g o + 3m (B A D 4 F 4 ) m k - j m 2 n C D2 {tt ) 0, m Si , m--, m s o) j 



+ 
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jj=l 




1=1 



Figure 13: Vertex correction with chargino-neutralino-w loop. 



with 



B 4 = 
C A = 
D 4 



3 cos 6 W u 



■ u N 



Zg (\zf sm9 w + Z 2 J* coaOw) + Y* z£ +3)i Zf 



sin6> 



w 



e rvMi* rvll , V M v (M+3)i* v 2l \ n MJ 
■n Z U Z + + Y u^U Z + I ^ 

sin t>\y 

J rjMi* ry2l*riMJ 



F 4 = -Y d J Z^*Z z _ l *C 



(B.12a) 

5 7 (B.12b) 

(B.12c) 

(B.12d) 

(B.12e) 
(B.12f) 



Diagram 5: Box with u 
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dr 5 dr tree _ MxpA vtt(uv) 



dAVdoosfl,,, dE e dcos9 eP ~ Sn ~ V ^ e ^p y 4m n f/ 2 \/ 2 16tt 2 

(jm (A 5 D 5 F 5 G 5 ) m e ) £^(0, 0, 0, m fij , m^, m«o , m~-)6' 

+ 2 3m (B 5 D 5 E 5 G 5 ) m v m k -m~^ n D (0, 0, 0, m-- , m s ., m^, m s o ) 

+ (jm (B 5 D 5 F 5 G 5 ) m- K -m 2 n + 3m (B 5 C 5 E 5 G 5 ) m s om^ 
Ddi(0, 0, 0, m-- , , , m«o ) (B.13) 

+ 3m (B 5 C 5 F 5 G 5 ) m 3 n D D2 (0, 0, 0, m k - , , m Vi , m s o ) 
+ . . . 
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with 



A 5 = - YfZ^Z 21 - 



B R = - 



sin 9w 

^--[—^ Z u Z + +Y u Z v Z + jC 
D§ — — Y J Z M ^* z 2l *C MJ 

_ [ 2^/2e (I+3)j yln v ly4nylj\ cMI 



MJ 



G, = 



Z\j (^Z 1 ™ sm0 + Z 2 ™ cos 9 

-= Z^* (Ztf sin 9 - Z 2 £ cos 9) 

V2sin0cos0w 



V^sin^cos^ 



\rl rvAn* rv(I+3)j 

1 u Zj N 



Diagram 6: Box with d 



(B.14a) 
(B.14b) 

(B.14c) 
(B.14d) 
(B.14e) 

(B.14f) 
(B.14g) 



jj=i 





i=i 



,K=1 



-,L=1 



Figure 15: Box diagram with d in the loop. 
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with 



dE e d cos9 e v dE e dcos9 



dr tree _ fp e xp p \ m 2 w <Re(UV) 

Sr 



E e E p J Am n U 2 V 2 16tt 2 



am (A 6 DqF 6 G 6 ) m e )D lu> (Q, 0, 0, ra d ~. ,m e - i ,m-+,m s o)5' i1 ' 

+ 2 am (B & D & E & Gq) m p m k om- K +D (0, 0, 0, m f q,mg. 1 m hl m- K +) 

+ (am (B 6 D 6 F 6 G 6 ) m k om 2 n + am (B 6 CqE 6 G 6 ) m k +m 2 n ^ 
D m (0, 0, 0, m^o, m^, m gi , m-+) 

+ am (B 6 C 6 F e G 6 ) mlD D2 (0, 0, 0, m s o, m^, m gi , m-+) 
+ . . . 



(B.15) 



ZJU I T/"A ^Ot XV 
/, , jv "+* Y i ^N^L 

COS fc'v 



nv 



B 6 



v^2 sine cos 9w 
—e 



Zf* (Zf sin 9 + Zf cos 0) + YfZ^Z 



K y3l* y(K+3)i* 



(B.16a) 
(B.16b) 



6 



sin 6 1 cos 9w 

D i" d TV D 



y) | sin^ - Z^costf ) f YjZ%** )j Z% ] 5 MJ (B.lGc) 



-'(J+S)j ylZ* , V J y3l* yJj 1 rMJ 



3 COS 9 W 

E« = - Y M Z^*Zl n C IM 



F 6 
G 6 



u 

e 



sin^ D 



d D 



glj* gin* _|_ -^/^(^+3)j^2n* j (jl^l 
e —Liyln | v L 7 (i+3)i 7 2n 



sin^ 



ir 



i L 



(B.16d) 

(B.16e) 
(B.16f) 

(B.16g) 



C Loop Integrals 

In this Appendix we have summarized the definitions of the loop integrals which were used in 
the calculation of the Feynman diagrams in the previous Appendix. 

C.l 2— point Functions 

First we define the 2-point functions as in because most 3- and 4-point functions will be 
expressed in terms of these. 
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with 



B (q , mi, m 2 ) — A— da; In if 

'o 



1 f 1 

Bi(q 2 , mi, m 2 ) = — -A + / dxxlnii 

2 Jo 

i Z" 1 

3 2 i(q 2 ,m 1 ,m 2 ) = -A - / dxa; 2 lniJ 

= -A - / dxx(l -a;) In if 
6 7o 



if = [(1 — rr)m 2 + ijtij — x(l — x)g 2 — ie] 
A = - - *y E + ln47r 



(C.l) 
(C.2) 
(C.3) 

(C.4) 



(C.5) 
(C.6) 



where 7^ = — d ln fJ x ^ = 0,577216... denotes the Euler constant. Since the one-loop 

x=l 

corrections to D are finite, the terms oc A cancel in the combinations of B— functions that 
will appear below. 



C.2 3— point Functions 

Besides the standard integrals Cq and C^ u with vanishing external momenta, we need the 
following integrals which contain combinations of the Feynman parameters in the nominator. 



x + y 



pl pl—x 

C D i(0,0,mi,m 2 ,m 3 ) = dx dy- 

Jo Jo m((l-x-y) + m z 2 x + m\y 

1 



mo — mf, 



(B (0, m 2 , mi) - B (0, m 3 , mi) - B 1 (0, m 2 , mi) + B 1 (0, m 3 , mi)) (C.7) 



1 — x — y 



rl rl—X 

C D2 (0,0,mi,m 2 ,m 3 ) = dx dy- 

Jo Jo mf (1 - x - y) + m^x + m^y 

= —5 9 ( Bl (°> m 2, mi) - £1 (0, m 3 , mi)) (C.8) 

Z" 1 /' 1 ~ :E (1-^-y) 2 

C D3 (0,0,m 1 ,m 2 ,m 3 ) = dx d 2/^r; r~ 5 — ; 2~ 

Jo Jo mi(l-x-y)+mix + miy 

1 



m, — ms 



(5i (0, m 2 , mi) - 5i (0, m 3 , m x ) - 5 3 (0, m 2 , m x ) + £3 (0, m 3 , m x )) (C.9) 
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C.3 4— point Functions 

Here the standard integrals Dq and D^ v with vanishing external momenta are needed. In 
addition the following two integrals with Feynman parameters appear: 



D D1 (0, 0, 0, mi, m 2 , m 3 , m 4 ) 

-1 rl—x fl—x—y 



— dx dy dz 



x 



o Jo Jo [ml (1 — x — y — z) + m\x + m\y + mfz] 2 

'Si (0,m 3 ,m 2 ) - -Bi (0,m 4 ,m 2 )) 



4 — m^ Vmg — ml 



- 2^ 2 (5i (0, m 3 , m 2 ) - Ex (0, mi, m 2 )) ) (CIO) 
^2(0, 0, 0, mi, m 2 , m 3 , m 4 ) 

/•l z-l-s rl-x-y x 2 

— dx dy dz 2 

Jo io Jo [ml (1 — x — y — z) + m\x + m\y + m|z] 

= o f — 5— 2 ( 5 i (°> m 3, "^2) - B\ (0, m 4 , m 2 ) 

m| — m\ \m\ — m\ 

- B 3 (0, m 3 , m 2 ) + B 3 (0, m 4 , m 2 ) ) 
l _ ^ 2 (^i (0, m 3 , m 2 ) - 5i (0, mi, m 2 ) 

- 5 3 (0, m 3 , m 2 ) + 5 3 (0, mi, m 2 ) )) (C.ll) 



m\ — m| 
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